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EFFICIENT  EVALUATION  OF  POLYNOMIALS  AND  EXPONENTIALS 
OF  POLYNOMIALS  FOR  EQUI-SPACEO  ARGUMENTS 

INTRODUCTION 


The  evaluation  of  polynomials  at  equl-spaced  arguments  Is  a  recurring 
task  that  arises  In  many  applications.  When  a  k-th  order  polynomial  is 
written  In  nested  form,  its  evaluation  generally  requires  k  additions  and  k 
multiplications  pt  each  argument  of  Interest.  For  a  set  of  equi-spaced 
arguments, -we-w+H' demonstratei-that  the  multiplications  can  be  entirely 

r 

circumvented  (except  during  Initialization)  and  that  a  recursive  procedure 
employing  only  k  additions  per  stage  will  suffice  to  generate  the  sequence 
of  polynomial  values. 


For  an  exponential  of  a  polynomial,  an  even  greater  savings  is  possible; 
namely,  the  exponential  can  be  circumvented  (except  during  initialization), 
and  only  k  multiplications  per  stage  are  required  in  a  recursive  procedure. 
Memory  storage  is  also  kept  at  a  minimum. 
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EVALUATION  OF  POLYNOMIAL 

The  procedure  Is  best  introduced  by  way  of  example.  Suppose  we  want  to 
evaluate  third-order  polynomial 

Pj(x)  *  a  +  Bx  +  +  ux3  (1) 

at  the  set  of  equi -spaced  arguments 

xn  *  Xo  +  nA  for  n  -  0,  1 ,  2,  . . .  .  (2) 

That  is,  we  are  interested  in  the  values 

2  3 

03(n)  =  P3(xn)  *  a  +  bn  +  cn  +  dn  for  n  «  0,  1,  2 .  (3) 

where 

a  -  a  +  Bx0  ♦  yxJ;  +  yx3  , 
b  =  A(0  +  2yXQ  +  3yX2)  , 

C  =  A2(y  +  3yXQ)  , 

d  -  A3y  .  (4) 

To  this  aim,  define  difference 

02(n)  =  Q3(n)  -  Q3(n-1)  =  b  -  c  ♦  d  +(2c  -  3d)n  4-  3dn2  .  (5) 

Also  define 

Q1 (n)  =  Q2(n)  -  Q2( n-1 )  «  2c  -  5d  *  6dn  ,  (6) 
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and  observe  that 


Q^n)  -  01  (n-1 )  =  6d  . 

These  last  three  recursions  together  read  as 


(7) 


0-j  ( n ) 

02(n) 

Q3(n) 


0-j  (n-1 )  +  6d 
Q2(n-1)  +  01 (n)  * 
Q3( n-1 )  +  Q2(n) 


for  n 


2. 


(8) 


and  require  only  3  additions  for  each  n,  with  no  multiplications  whatsoever. 
The  starting  values  for  recursion  (8)  follow  Immediately  from  (6),  (5),  and 
(3),  respectively: 


Q] (0)  =  2c  -  6d  , 

Q2(0)  =  b  -  c  +  d,  (9) 

03(0)  =  a  . 

Extension  to  a  k-th  order  polynomial  Is  obvious,  and  requires  k 
additions  per  stage,  with  no  multiplications. 
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EVALUATION  OF  EXPONENTIAL  OF  POLYNOMIAL 

Suppose  we  want  to  evaluate  the  exponential  of  a  third-order  polynomial, 
namely 

2  3 

P3(x)  *  exp[a  +  0X  +  yx  +  yX  ]  (10) 

at  the  arguments  listed  In  (2).  That  Is,  we  want  the  values 

2  3 

03(n)  =  P3( xn)  =  exp[a  +  bn  cn  +  dn  ]  for  n  =  0,  1,  2  ...  ,  (11) 

where  a,b,c,d  are  given  in  (4). 

To  accomplish  this  goal,  define  ratio 
Q2(n)  =  03(n)/Q3(n-l )  *  exp[b  -  c  +  d  +(2c  -  3d)n  ♦  3dn2]  .  (12) 

Also  define 

Q1 (n)  =  Q2(n)/Q2(n  -  1)  =  exp[2c  -  6d  +  6dn)  ,  (13) 

and  observe  that 

01 (n)/0] (n  -  1)  =  exp[6d ]  .  (14) 

These  last  three  recursions  together  read  as 
Q^n)  =  0] (n  -  1)  exp[6d]  "'j 

Q2(n)  =  02(n  -  1)  Q^n)  for  n  =  1,  2 .  (15) 

03(n)  =  Q3(n  -  1)  Q2(n) 
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and  require  only  3  multiplications  for  each  n,  with  no  additions  or 
exponentiations.  The  starting  values  for  recursion  (15)  follow  immediately 
from  (13),  (12),  and  (11),  respectively: 

Q-j  (0)  =  exp[2c  -  6d]  , 

Q2(0)  =  exp[b  -  c  +  d]  ,  (16) 

Q3(0)  =  exp[a]  . 

Initialization  requires  the  evaluation  of  four  exponentials. 

Extension  to  an  exponential  of  a  k-th  order  polynomial  is  obvious,  and 
requires  k  multiplications  per  stage.  Initialization  requires  the 
evaluation  of  k  +  1  exponentials. 
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CONCLUSION 


All  the  results  above  apply  to  complex  coefficients  a,b,c,d  as  well  as 
complex  arguments  xQ,a.  Only  integer  n  needs  to  be  real.  However,  since 
a  complex  multiplication  involves  4  real  multiplications  and  2  real 
additions,  the  time  of  execution  will  naturally  be  larger. 

The  applicability  of  the  above  results  to  linear  frequency-modulation 
with  Gaussian  amplitude-modulation  follows  readily,  by  restricting  the  order 
of  the  polynomial  in  (10)  and  (11)  to  k  =  2;  i.e.,  set  y  =  d  =  0.  This 
particular  case  has  been  treated  in  [1];  in  particular,  the  accuracy  of  the 
procedure  has  been  investigated  and  found  adequate  for  most  appl ications . 

The  evaluation  of  cosines  or  sines  of  real  polynomials  can  be  achieved  by 
setting  a,b,c,d  In  (11)  to  purely  Imaginary  values. 

When  k  equals  2  and  a,b,c  are  complex,  the  quantities  Q^n)  and 
Q2 ( n )  are  complex.  Since  a  complex  multiplication  involves  4  real 
multiplications  and  2  real  additions,  the  number  of  operations  per  stage  to 
generate  Q^n)  and  02(n)  Is  8  real  multiplications  and  4  real 
additions.  As  an  example,  if  a,b,c  are  purely  imaginary,  a  =  i a ' ,  b  =  ib', 
c  =  ic ' ,  then 


Q.,{n)  =  exp [ i ( a '  4-  b'n  4-  c'n^)]  = 

2  2 

=  cos[a'  4-  b'n  +  c'n  ]  4-  i  sin[a'  4-  b'n  4-  c'n  ]  ,  (17) 

meaning  that  the  cosine  and  sine  are  capable  of  simultaneous  generation  at 
each  stage.  (Attempts  to  generate  the  cosine  alone,  with  a  lesser  number  of 
operations  per  stage,  have  not  been  successful.) 
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There  1$  no  need  to  set  aside  storage  arrays  for  the  recursive 
quantities  in  (8)  or  (15),  if  these  numbers  are  used  on  the  fly  as  they  are 
generated.  Then  the  computer  coding  for  (8)  is  simply  Q1  =  Q1  +  06, 

Q2  *  Q2  +  Q1 ,  Q3  =  Q3  +  Q2,  while  that  for  (15)  is  simply  Q1  =  Q1*E6,  Q2  = 
02*Q1 ,  Q3  -  Q3*Q2,  in  the  order  listed.  Generally,  storage  of  only  k 
temporary  variables  is  required  for  a  k-th  order  polynomial.  On  the  other 
hand,  if  03(n)  must  be  stored  for  later  use,  the  only  change  in  the  coding 
above  is  to  replace  the  Q3  lines  by  Q3(N)  =  Q3(N  -  l)  +  02  and  Q3(N)  = 

03 (N  -  1)*Q2,  respecti vely;  there  is  no  need  to  store  Q1  or  Q2  in  arrays. 

For  general  values  of  k,  if  the  product  of  a  rational  function  with  an 
exponential  of  a  polynomial  must  be  calculated,  it  can  be  broken  down  into 
the  evaluation  of  two  polynomials  and  one  exponential,  as  indicated  above. 
Then  one  additional  multiplication  and  division  realizes  the  desired 
combination.  Extensions  to  sums  and  products  of  such  combinations  are 
obvious.  The  orders,  k,  of  the  numerator  and  denominator  polynomials,  as 
well  as  the  polynomial  inside  the  exponential,  need  not  be  equal,  but  are 
completely  arbitrary. 
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